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DYNAMIC BEHAVIOR OF AN UNSTEADY TURBULENT BOUNDARY UYER 


P. C. Parlkh, W. C. Reynold*, R. Jayaraann, and L. W. Carr 
Dapartaant of Maehanleal Entlnearlng, Stanford Univaraity 
Stanford, California 9A30S 


Suaaary 

Thia paper raporta axparlaanta on an unataady turbulent boundary 
layer. The upatraaa portion of the flow ia ataady (in the naan). In 
the downatreaa region, the boundary layer aaea a linearly decrcaaing 
free-*atreaa velocity. Thia velocity gradient oacillate* in tine, at 
frequenciea ranging froa aero to approzinately the burating frequency. 
Conaiderable detail ia reported for a lo«^aaplitude caaa, and pral^inary 
results are given for a higher aaplitude sufficient to produce aoae re- 
verse flow. For the snail anplituda, the nean velocity and naan turbu- 
lence intensity profiles are unaffected by the oscillations. The 
aaplitude of the periodic velocity conpooent, although as much as 70Z 
greater than that in the free streaa for very low frequencies, becones 
equal to that in the free streaa at higher frequencies. At high frequen- 
cies, both the boundary layer thickness and tl.e Reynolds stress distribu- 
tion across the boundary layer becoaa frozen. The behavior at larger 
aaplitude is quite similar. Most iaportantly, at sufficiently high fre- 
quencies the boundary layer thickness mains frozen at its mean value 
over the oscillation cycle, even though flow reverses near the wall during 
a part of the cycle. 


Introduction 

The objectives of the Stanford Unsteady Ihrbulent Boundary Layer 
Prograa are: to develop a fundamental understanding of such flows, to 
provide a definitive data base which can be used to guide turbulence aodcl 
developaent, and to provide test cases which can be used by cc»putors for 
coaparison with predictions. 

Due to space liaitations, work of other investigators will not be 
suaaarlzed here, except to note that all the previous experlaents are 
characterized by unsteady flow at the inlet to the unsteady region. For a 
comparison of the present experimental parameter range with those of other 
investigations, see Reference 1. The distinctive feature of the present 
experiments is that the boundary layer at the Inlet to the unsteady region 
is a standard, steady , flat-pXate turbulent boundary layer. It is then 
subjected to controlled oscillations of the free stream. This feature ia 
especially important from the point of view of a coaputor, who needs pre- 
cise specification of boundary conditions for computation of the flow. 
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Pf-StfMi Boundary Coodition of th> Pr«««at B»pTl««ot 

Th« dMircd fr««-strMB v*loelty u«(x,t) in tlM waur t'umal built 
for this work la ahovn In Fig. 1. u« raaalaa ataady and ualfore for tba 
flrat two aatara of boundary layar dayalofUMat* It than -daeroaaoa lln- 
aarly In tba taat aacclon; tba nagnltudo of tba voloelty gr.adlont varlaa 
alnuaoidally fron taro to a naxiaun valua during tba oaclllatlon eyela. 
Tba naan frar*atroan voloelty dlatrlbutlon in tba taat aaetlon la tbua 
llnaarly dacraaaing and corraaponda to tba dlatrlbutlon at tba eyela pbaaa 
angla of 90* , whlla tha anplltuda of Inpoaad fraa-atraaa oaelllatlona 
growa llnaarly in the atraanwlaa dlraetlon, atartlng at laro at tha 
aotranca to a aazinun valua of a^ at tba axit. Honea, 

Ujx.t) - * < % 

• “-,o ■ - • *0 < * < % 

The iaportant paranatara of thla problan ara tha aaplituda pa.vaaatar 

a ■ q and tha fraquaney paranater: 0^ ■ ffi^/u. Hara f ■ 

u/(2«) and 6^ is the thieknasa of tha boundary la^ar at tha inlet to 
the unsteady region. In tha prasant azpariaanta: 

u, o • 0*^3 a/s, 6 ^ - 0.05 a, 0 <f < 2 hz, 0 < o < 0.25, 0 < < 0.14 

It should be aantioned that tha valua of tba fraquaney paraaatar 0^ 
at the ao-callad "bursting fraquaney" in turbulent boundary layers la 
about 0.2 [2]. Thus the laposad oscillation frequancles used in tha 
present expariaants cover tha range fron quasi-steady (f ■ 0) to values 
approaching the bursting frequency. Tha rasi'lts reported hare arc for two 
non-diaenslonal aaplltudes, a “ 0.05 and 0.25 (ncmlnally). The latter 
is sufficient to cause reverse flow in a turbulent boundary layer at the 
end of the test section during a part of the csclllatlon cycle. 

Experiaental Facility 

Figure 2 is a acheaatlc of the facility. Tba 16:1 nozzle contraction 
is followed by a 2 a long dcvelopncnt section, where the test boundary 
layer is grown on the top wall. A constant head and a constant flow re- 
sistance provide a constant flow . The free-streaa velocity in the devel- 
opnent section is naintalned unifora along x bv bleed fron the bottom 
wall. 

The linear decrease in free-streaa velocity in the test section Is 
accenplished by uniformly bleeding off sooc flow through the botton wall 
in the test section. The rouinder of the flow exits downstreaa. Each of 





chpM two flows oxita cho tu&nol through slots In sn oscillstlng plsts. 
Ths daslgn snsurss that, rsgsrdlsss of tbs position of tbs oscillstlng 
plsts, .ths total flow arsa of ths slots raaains ths saaa. Iho slots sro 
ths eontroUlng raslstaneo of ths antlra fluid circuit, banco ths constant 
flow. By sinusoidally oscillating tho plats, a linaarly docrsaslng peri- 
odic frss-stroaa distribution is astahllshod la tbs tost ssctlon, whlla 
ths upstrsaa flow in ths dswolopaont ssctlon roaalns steady. 

Msasursaant sod Data-Procssslng TOchniquas 

Pitot tubas sra used for naan waloclty Masuranants in steady flow 
raglons. Unsteady velocity aaasureuants use a single-channel, forward- 
scatter, Bragg-shifted DlSi laser anaaoaeter la ths tracking nods. 

Following Ikissaln and Beynolds [3], the instantaneous velocity signal 
from an unsteady turbulent flow aay be decouposed Into three parts: 

u-u+u+u' (1) 

where u is the aean, u is the tine-dependent, organised (detemlnlstlc) 
coapooent, and u* la the randon fluctuation, u is deteralned by long- 
tine averaging of u. Here u is of a periodic nature and aay be deter- 
nined by first phase-averaging the instantaneous velocity signal and then 
subtracting out the aean. Thus , 

S - < u > - u (2) 

Here < u >, the phase average velocity, is deteralned by averaging over 
an enaenble of saaples taken at a fixed phase in the inposed oscillation. 
In the present expetiaents, with haraonlc oscillation of the free streaa, 
the response at points within the boiuidary layer is alaost sinusoidal, 
with higher haraonica contributing less than 5Z. Hence, u aay also be 
extracted froa the instantaneous signal u by cross-correlation with a 
sine wave In phase with the oscillation. A digital correlator (HP 3721A) 
was used to deteralne cross-correlations leading to the u data reported 
here. Currently a DEC MINC-IL laboratory ainicoaputer systea is used for 
autoaatlc data acquisition and processing, allowing the deteminatlon of 
phase averages of u and u* . 

The aeasureaents reported here were taken at a fixed streonwlse 
location near the end of the test section at x - x^ ■ 0.568 a. 
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Boh>vior at Saall A«plitud« of Lipo««d OtclllatioM 

The MM velocity profllct MMurod irith the eoelllatliit plat* Ic 
fixed poaltioDa 6 • 0, 90*, 180* ere fit Vy daahad eurvea in Fig* 3* 
Theaa repraaeat phaae*averaged prefilea at aero frequeacy, i.e*, quaal* 
ateady profilaa. At thia aaplltude (a ■ O.OS), tbe reapoaae of the 
bouadary layer la alMat llaear , ao . that the profile oorreapoadlag to 0 ■ 
90* liea nearly aldway betweea the 6 ■ 0 and 180* profilaa. The 90* 
profile repreaenta the aeM profile for quaal-atoady oaeillatlona . The 
difference betweM the 0 and 90* profilaa at a fixed y-locatlon 
repreaenta the aaplltude of quaal-ataady oaeillatlona at that location in 
the- boundary layer. Note that the quaai-ateady.aaplltudea in the boundary 
layer are larger than the free-atreaa aaplltude. 

The aean velocity profl.lea aeaaured under oacillatory condltiona at 
0.5 hz and 2.0 hz are ahovn aa data polnta in Fig. 3. Note that the 
aeM velocity profllea at varloua frequanclea are identical with the pro- 
file Manured under atationary condition with pulaer angle aet at 0 ■■ 
90*. It My be concluded that the Man velocity profile (at a fixed aa- 
plitude a " 0.05) ia Independent of tbe lapoaed oacillatlon frequency in 
the entire range 0 ^ f ^ 2 hz. The aaM behavior peralata all the vay up 
to the vail. 

Thia behavior of the Man velocity profile My be explained by an 
exMination of the governing equationa. Uae of (1) in the aoaentia equa- 
tion and tlae-aver aging ylelda 

Equation (3) My be recognized aa tha equation governing an ordinary tur- 
bulent boundary layer, except for tbe addition of the tera uv', which 
repreaenta Reynolda atreaaea arialng free the organized oaeillatlona. 

nie tiae-eean prcaaure gradient 3p/0x My be ahovn to be indepen- 
dent of the lapoaed oacillatlon frequency and the aeM aa that obtained 
for f - 0 at 6 « 90*. Therefore, the Man velocity field will be 
frequency-dependent if and only if one or both of the following happen: 

• The diatributlon of Reynolda atreaa u'v* ia altered under oacllla- 
Cory condltiona and ia dependent on the frequency of lapoaed oacll- 
latlona. 

• The Reynolda atreaa uv a rialn g froa organized fluctuatlona becoaea 
aignlf leant coapared with u*v' . 


He shall DOW argus that nalthar sf tha abova ra^ulraMnts la aac. Figuta 
4 sbova tha aaasurad diatributiea of uadar stationary condition 

with ths yulsar at 9 • 90* ( tha naan position) as wall as thosa naa- 

surad tuidar oacillatory conditions at fraquancias up to 2 hs. Nota that 
is indapandant of tha iaposad oseiUatlen fr^uaney and, furthar, 
that it is the Sana as that naasurod at f ■ 0 and 0 ■ 90*. Ha baliava 
that tha saaia would ba trun for which at praaant wa cannot naa- 

sura. Figura S gi^aa a eonparison batwoan naasurad waluas of uv at 2 hs 
with data on u’V’ obtainad by Andaraon {4) in a staady advaraa prassura 
gradiant boundary layar at conparabla conditions. Tha prasant data on 
uw wars obtainad by -saparata L&4 nanausanants of u and a and thair 
raspactiwa phasas. It nay ba aaan that tha contribution of uv to total 
ftaynolds stress is insignificant oaar alaost tta, satire boundary layar. 
Banca, is indapandant of frequency and us is nagligibla, and so 

tha naan velocity profile is also indapandant of frequency and is tha sane 
as that found at f ■ 0 with 6 ■ 90*. 

Tha behavior of the periodic conponant u will next ba exaninad. Ha 
denote 

^ ■ •l<y) coalnt + ♦(y)] (4) 

Tha profiles of aaplitudas Sj^ aaasurad in tha boundary layar and noraal’* 
ixad by the free-streaa aaplitude a^ ,, are shown in Fig. 6. Tha profile 
for quasi-steady (f ■ 0) oscillations was datarainad, as explained ear- 
lier, froa tha aaan velocity profiles naasurad at f ■ 0 with 6*0, 
90*, and 180* (sea Figs. 3(ay , (b)}. Note that, during quasi-steady 
oscillations, tha aaplitude in tha boundary layar exceeds tha frea-streaa 
aaplitude by as auch as 70Z. It nay be aantioned that data for f ■ 0. 1 
hx, not shown on Fig. 6, do indeed coaa vary close to the quasi-steady 
behavior. 

As Che frequency is increased, the aaplitude within the boundary 
layer is attenuated. The aaplitude appears to drop as f is increased 
and then rise again. At high frequencies, the aaplitude in nost of the 
boundary layer is the saae as in the free-streaa; near the wall the aapli- 
tude of the periodic coaponent rapidly drops to sero. 

The phase differences between the boundary layer oscillations and 
free-atreaa oscilli clone are shown in Fig. 7. For f * 0 there is no 
phase difference. The largest phase lags in the outer region of Che 
boundary layer were observed at f * 0.25 hz. The effect of increasing 
the frequency is to reduce the phase lag in the outer region, bu*^ 
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Introduce Inrt* plmse loodo In thn rnsion unrp elOM to the mU. 
Clonrly, the MToptotie bobaulor of tho outer rotiw for hith frequeneiee 
le once egeln e lero pheee leg with reepeet to free^etreea oeellletlone • 
ee in tla^queel^eteady'ceeo* 

At high freqoenclee, the nonblnetlon of the eeynptotle beheriore of 
e|/e^ . end d In the outer reglM together -oith-the^ect thet the Men 
velocity profile le uneffected by Inpoeed oeellletiene , bee the effeet..of 
freeslng the boundery leyer thlckneee. Thia -le ehoun la Fig. 8, where the 
phaee-avereged boundary layer thickaeea < 8^99 > la plotted aa a func- 
tion of ^he -cycle, pbaae angle, for. eeweral. freqoeneloa.. The quaalrateady 
behavior of - < > la quite -obvlouat at 6 • 0, the . boundary, layer 

in the teat aeetion concinuea to develop -under a aero preaaure gradient 
and la the thlnneat at thla point in the entire cycle. Aa the phaae angle 
la Incraaaed, preaaure gradienta of increasing adveraity are inpoaed on 
the boundary layer, cauaing it to thicken. The naslnun thlckneaa ia at- 
tained at 6 ■ 180* under the nazlnue adverse pressure gradient. Hence, 
at f ■ 0, 6 oacillatea 180* out of phase with u.. 

Under oscillatory condltlona at f - 0.25, 0.5, and 2.0 ha, two 

things happen: a significant phase lag develops fro* quasi-steady behavior 
and the anplitude attenuates with increasing frequency. For the f ■ 2.0 
hs case, the verletion over the-con p lete-cycln— is less than IZ and the 
boundary layer thickness In practically frozen during the oscillation 
cycle. 

It nay be abown by a sinple argMsnt based on a nixing length nodal 
of boundary layer turbulence that the freezing of the boundary layer 
thickness at high frequencies is also acconpanied by freezing of the Bay- 
oolds stress over the oscillation cycle, lb prove this, we hypothesize 
that the phase-averaged Beynolds stress distribution My be related to the 
phase- averaged velocity profile 
ary layer, i.e., 

- < u*v' > - c 

a 

Now, in the outer region of the 
be nodeled as 

i 

where 1 is nearly a constant. 

< u > • u u 


In the seM Mnner ae for a steady bound- 


a < u > 

“■57“ 


C - 1" 




(5) 


boundary layer, the nixing length 1 nay 


Now, 


( 6 ) 


■ 11+ Sj^(y) cos[ut + i|»(y)] 


( 7 ) 


-h- 





Hovctm, 1b tha bl(lrfr«9ttBBe7 11 b1c» 

•^(y) “ " coBBt ; ♦(y) • 0 BBd < « > • T ■ cooBt. 

Thar* for* 

a < tt > _ iu 

“TT“ ^ 

CoabiniBg cba abetrai ooa floda 



( 8 ) 

(9) 


( 10 ) 


l.a.| tha phaaa-aTaragad laTiiolda atraaa la tba outar ragloa alao baeoaaa 
froaaa at •> u'v*. 

Expariaantal avldaaea of thia atraaa- fraaalng bahavlor vaa obtainad 

2 

bp BaaBurauata of phaaa-avaragad aomal turbulant atraaa < u* >• Tba 
quaal-ataady (f - 0) profllaa of < u'^ > ara ahova In Fig. 9 for thraa 
phaaa anglaa 6-0*, 90*, and 180** Not a that tha dlatrlbutlon for 90* 
Ilaa naarly oddway batvaan thoaa for 0* and 180*. TIm d latr lbutlon of 
< u'^ > for 90* la tha aaaa aa tha dlatrlbutlon of u* , aa aaan 

earllar. Tharafora, tba dlffaranea batwaan tha 0* and 90* curvaa In 
Fig. 9 rapraaanta tha aaplltuda of quaal-ataady oaclllatlona of < u'* > 
at any point in tha boundary layar. Thla aaplltuda waa dataralnad graptr 
Ically froa Fig. 9 and la plottad la Fig. 10 for tha eaaa of f • 0. 
Undar oaelllatory eendltlena, tha aaplltuda of tba noraal atraaa oaell- 
latlona In tha boundary layar attanuataa aa tha fraquancy of lapoaad oa- 
clllatlona la Incraaaad froa f - 0. At f ■ 2.0 ha, tha aaplltuda of 
atraaa oaclllatlona aeroaa tha boundary layar la alaoat aaro oaar tha 
outar raglon, aa aaan la Fig. 10, l.a., tha atraaa la alaoat froaan ovar 
tha oacillatlon cycla. 


Bahavlor Undar Larga Aaplltudaa of tepoaad Oaclllatlona 

Ha now diacuaa tha caaa of a - 0.25. All data raportad for thia 
caaa ara prallalnary and aubjact to rawlalon. They ara Included hare 
bacauaa of thalr apaclal Intaraat to thia naatlng. Alao, bacauae of 
apparatua pacullarltlaa , o warlaa aoaawhat with f In tStla caaa , 
banco 0.25 la only a noalnal walua. 

The bahavlor la qualltatlvaly alallar to tha o “ 0.05 caaa. Tha 
naan velocity profllaa for f - 0, 0.25, 0.5, and 2.0 Ha are ahown in 
Fig. 11. Nota that tha profllaa ara Identical for tha caaaa of f ■ 0, 
0.5, and 2.0 hz. For tha caaa of f - 0.25, however, there la a 
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> > 4g >H4c« itr ^ tta MMr part 'vf tti fUs da«l- 

AtlM r««tilt« fr«i «xe«Mlv<i thick— lag of the b—c dory layer dorlag * 
part of tha eaelllatioa eyela are— d the phace aagla ef 100*. The hi— k- 
ate eff—t ef- — -'am— lealy thick be— dary layer aa— — iae tease la 
the local free-stre— eoleelty la tha t— t a— t l e a* Thareferat the 

destrod ll—arly d— r— slat free*stre— eeloelty dlstributi— la — t 
aehiatad eter a part ef tha cycle. At -hither fraq— 1— •- theagh. tha 
b— adary. layer thick— as orer tha eatlre -eaciUatlea cycle datlates eery 
little fr— its aa— tal— . cerreapeadiat to the 9 ■ 90*, t • 0 coa> 
dltlon. 

llM-behatlor of tha a— lltuda ratio aad pha— difference with reapect 
to' free -etre— , as aheim ia yits. 12 aad 13, is quits slailar to that for 
the loaar*— plltude ca—. At high fraq— acy, tha orarahoot ia tha — pli- 
tuda ratio disappears and phase angles over —sc of tha boundary layer 
approach saro. Vary close Co Che vail, chare is a Candaney to develop 
phase leads. 

The phaae>averaged velocity profiles for f * 2.0 hs are shovn in 
Fig. 11. HoCe Chat at 8 • 180* Chare is a aaall region of reversed flow 
close to Che vail. Despite this flow reversal, the boundary layer chick- 
nesa re— ins close to its —an value, — seen in Fig. IS. This behavior 
is in conefast to that of a steady boundary layer, where excessive thick- 
ening of Che boundary layer occurs as flew reversal is approached. At 
low frequency (f " 0.2S ha), the chick— as oscillates as «uch as ± 40Z 
about Che —an val— ; however, at f ■ 2.0 hs this variation is only 
about t SZ. 

Conclusions 

The conclusio— fr— our ezperlaents to date — y be sun— rlsed as 
follows: 

1. The —an velocity profile in the boundary layer is u— ffected by 
laposed frec-stre— oscillatlo— in the range of iTequencies en- 
ployed, and it is the sa— as the o— — — ured with a free-strean 
velocity distribution held steady at its mean value. 

2. This behavior of the —an velocity field is a co— equence of two 

observations: (a) the tlae<*averaged Baynolds stress distribution 

across the boundary layer is u— ffected by the inposed oscillatlo— 
and is indeed the sa— as the o— a— sured with the free-streaa vel- 
ocl':y distribution held steady at the a— n value; and (b) the Rcy- 
— Ids stresses arising fr— the organised velocity fluct— tio— undur 
iapoaed oscillatory conditio— are — gligible c— pared to the Rs]r- 
— ids stresses d— to the rend— fluct— tio— . 
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3. Th« Miplltud* of eh* p«riodie ea^Mmt la tht boaadar^ l&ft aadar 
quaal-ataady eaclUatloaa (f ^ 0) la aa aaeh aa 70S largar than tha 
lapaaad fraa-atraaa aaipllttada. BMiavar, ac hlglMir fragnaaelaa tha 
pa^ aaplltttda la tha bovadarp lajar la rapldlp attaauatad towrd aa 
aajaptetie bahavlor vhara ai^ltadaa la tha oatar T^iea of tha 
hooadary lafar h ac oaa tha aaao-aa tha fraa-atraaa aaplltiido, dropplag . 
off to goto la tha aaac-wall ragloa. 

4. Quaai-ataadp houadarp lapar valoeltp raapoaaa la la phaaa with tha 
i^aad .fraa-atroaa oacillatioaa*; Aa tha fraqaaaep la iaeraaaod, 
phaaa laga bagla to dawalop la rtha . oatar ragioa of tha bouadarp 
lapar. Iho aagaltuda of ttia ^Moa lag roachoa a aaalaua aad thaa 
daoraaaaa with iac;aaalng Jrogaaaep uatll -as aapaptotle Halt la 
raachad vhara tha oatar r^loa oaeo ^ala raapoada la phaaa with tha 
fraa atraaa. Hoar tha Mil, howawcr, largo load aaglaa ara praaaat 
at hlghar oaeillatloa fraguaaelaa* 

5. A eoeaaquaaea of <3) aad (A) abowa la that tha bouodarp lapar thick- 
oaaa baeoaaa aaarlp froaaa owor tha otclUatloa epela at hlghar fra- 
qoaaclaa. Thla raaalaa trua awaa If flow rawaraal takas placa la tha 
naar-wall ragioa ovar a part of tha eaciUatloB epela , aa la tha 
larga-aoplltuda eaaa. 

6. A eoasaqaaaea of (3}. (4)» and (5) abowa la that tha lapaolda atraaa 
dlatributioa la tha outar ragioa of tha bouadarp lapar also baeoaaa 
froaaa orar tha oaeillatloa epela at hlghar fraquaacias. 
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